Oscillatory screening of the dc electric field in the Si-Si02 
multiple quantum wells probed by second-harmonic generation 



V. V. Savkin, A. A. Fedyanin, A. N. Rubtsov, and O. A. Aktsipetrov 
Department of Physics, Moscow State University, 119899 Moscow, Russia 

F. A. Pudonin 

Lebedev Physical Institute of RAS, Leninsky pr. 53, 117924 Moscow, Russia 

(February 1, 2008) 

Abstract 

DC-electric field, being screened in 3D semiconductors, normally decays 
monotonically in space. Experimental studies of the DC electric field screen- 
ing in Si-Si02 multiple quantum wells by electric field induced optical second- 
harmonic generation show a non-monotonic, oscillatory-like decay. The model 
of electrons localized inside quantum wells, with the first subband occupied, 
allows a description of the phenomenon. Interwell Coulomb interaction, a 
finite value of the electron charge and strong effective-mass anisotropy result 
in a crucial difference from 3D Fermi liquid. 
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Properties of the degenerate electron subsystem of a solid are basically determined by an 
interplay between Coulomb interaction and quantum effects, which is revealed in correlation 
and exchange interaction. The corresponding typical energy scales are given by character- 
istic magnitude of the Coulomb interaction ec oc e^rr^l'^ and the Fermi energy for 
noninteracting electrons Sp (x fi^m^^n~'^^^ ^ where n is the electron density, G denotes the 
lattice dielectric constant, mo, e and Ti are fundamental constants. \i Spjec >> 10^^ 10~^ 
the system has the properties of a normal Fermi liquid [IH, whereas in the opposite case a 
dielectric bandgap forms 0. An Anderson-Mott transition occurs at the specific value of 
ep/ec i 

Significant attention has been paid, over recent decades, to spatially confined 
semiconductor-based systems. Advanced modern technologies allow the fabrication of struc- 
tures with a predetermined morphology, and therefore control over their electronic proper- 
ties. Amongst 2D semiconductor-based systems (superlattices and multiple quantum wells) 
two extreme cases should be pointed out. While epitaxial structures like GaAs/AlGaAs 
are characterized by relatively close bandgaps of the components, the bandgap difference in 
amorphous Si-Si02 based systems are much higher than the Si bandgap. 

The anisotropic dispersion law of electrons requires the consideration of two parame- 
ters, {eF/£c)\ \ and {eF/£c)±, which describe tangential and transversal electron motion, 
respectively. For the Si layer with a carrier density of 10^^ 10^^ cm~^, at the zero tem- 
perature and with m ^ rriQ we have {eF/£c)\ \ ~ 10^^ 10^^, i.e. lateral properties of the 
system can be described in terms of a Fermi-liquid. The parameter {sp/ec)^ is governed 
by the transparency of potential barriers in the superstructure. For high potential barriers 
the parameter {ef/sc)^ is very small, and the interlayer electron motion is suppressed by 
the Coulomb repulsion. Electrons are localized within the layers and the structure can be 
treated as a set of independent quantum wells (QW), i.e. as multiple QW's (MQW). The 
Mott-like transition is expected at a certain transparency of the barrier. 

Properties of the electron liquid reveal themselves in the response of the system to exter- 
nal fields. For example, Friedel oscillations appear in the screening cloud in the degenerate 
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Fermi-liquid [Q. The period of these oscillations is governed by the typical de-Broiglie wave- 
length in the system. OD systems weakly interacting with their environment show a Coulomb 
blockade phenomenon, which is due to a finite value of the elementary charge |p. Analysis 
of four coupled quantum dots shows that properties of an array of interacting quantum 
dots depend on the interdot tunnel coupling. For weak coupling this system behaves as a 
set of independent dots, at intermediate coupling the Mott insulator forms (the so-called 
collective Coulomb blockade) and for strong coupling the Coulomb blockade is completely 
destroyed. It is worth noting that for the collective Coulomb blockade the screening has 
an essentially collective nature, but it still shows non-monotonic features. Particularly, the 
carrier density in neighboring quantum dots is anti- correlated. 

On the other hand, the DC-electric-field screening in bulk semiconductors shows mono- 
tonic spatial distribution of the voltage drop across the space charge region . We conclude 
that nonmonotonic peculiarities in the DC-screening indicate that the electron subsystem 
essentially differs from the bulk 3D case. 

Optical second-harmonic generation (SHG) is a versatile optical probe of surfaces and 
interfaces [Q. This technique is applicable to all interfaces accessible by light, including 
buried interfaces 0, and has inherently high spatial |TD[, spectral [|rT[] and temporal 



resolutions. The sensitivity of SHG arises from the strong symmetric selection rules for the 
second-order susceptibility. Namely, the second-order susceptibility vanishes in the dipole 
approximation for media with inversion symmetry, and the quadrupole term becomes preva- 
lent. In the presence of the DC- electric field the dipole contribution into the second-order 
polarization P^"^ arises ||13||: 



where E'^ and E^'" are the fundamental and DC fields, respectively; k'^ is the fundamental 
wavevector; nonlinear susceptibilities x^"^^^ ^ind x^^^^ ^ire responsible for the 2lu = cu + uj and 
2uj = 0+uj+uj processes, respectively. Two terms in (|ip are comparable at E^'-' 10^ V/cm. 
This sensitivity to the DC-fields has been experimentally demonstrated for semiconductor- 



electrolyte and semiconductor-dielectric interfaces . 

In this paper, the DC electric field screening in Si-Si02 MQW is studied by DC-field 
induced SHG. The magnitude of parameter {eF/£c)x_ smaller than 10~^ for the structures 
studied. The experimental bias dependence of the SHG intensity for MQW reveals the essen- 
tially non-monotonic, oscillatory-like behavior. This is attributed to the essential quantum 
effects in the highly anisotropic layered system with a Coulomb interaction between layers. 

MQW structures were fabricated by RF magnetron sputtering on a vicinal n-Si(lOO) 
wafer [Q. MQW structures consist of 40 double amorphous layers of Si-Si02. Two MQW 
samples were studied with Si02 layers 3 and 5 nm thick. The Si layers were 1.1 nm thick 
for both samples. The perfect periodic layered structure of the MQW is studied by X- 
ray diffraction and on-line Auger-electron spectroscopy. The amorphous structure of the Si 



layers is confirmed by Raman spectroscopy |14 



The output of a Q-switched YAG:Nd^^ laser at a wavelength of 1064-nm was used as 
fundamental radiation. A pulse duration of 15 ns and an intensity of 1-10 MW/cm^ avoid 
laser damage of the sample and other undesired photo-induced effects. The experimental 
setup is described in detail elsewhere |T5|. The transparency of the Si substrate at the 
fundamental wavelength allows the study of SHG in transmission. The advantage of this 
geometry is that only a few components of the nonlinear susceptibility tensor contribute to 
the signal, which simplifies the analysis. 

DC-electric-field is applied to the MQW by imposing a bias voltage between the Al back- 
side electrode and the In-Ga ring gate electrode on the top of MQW. The free carriers easily 
move in the in-plane direction, while the QW's are well isolated from each other, therefore 
the voltage distribution inside the structure does not depend on the in-plane coordinates, 
although the cap electrode is not planar. 

The in-plane symmetry of the samples was analyzed by measuring the azimuthal SHG 
anisotropy. Figure 1 shows the dependence of the SHG intensity on the azimuthal angle tp. 
The two-fold symmetry of the dependence indicates the in-plane anisotropy of the sample 
attributed to the initial miscut of the substrate from the (100) direction. From the micro- 



scopic viewpoint, this means the "terrace" structure of the interface. We expect that these 
terraces are reflected in the morphology of several Si and Si02 layers, closest to the sub- 
strate, which possess the 2/m point-group symmetry (Figure 1). The deposition of further 
layers erodes the anisotropy. Therefore the symmetry of the rest of layers is expected to be 
oo/m. 

Since both 2/m and oo/m layers possess an inversion symmetry, P"^^ is given by Eq.(0). 
Symmetry of the layer governs the selection rules for the components of the fourth-rank 
tensors x'"^'^'' ^ind x^^^"* ■ Take the y-axis to be parallel to the terraces, and the 2;-axis normal 
to the surface. For the transmission geometry xzxx, xzyy, and yzxy = yzyx components 
of x*'^'^'' s-i^d x*-^"^-* contribute to the SHG intensity for the 2/m layers, and there is no 
contribution from oo/m layers. Consequently, the DC-electric- field effects in the buried 
Si-Si02 layers near the vicinal Si substrate are purely detected. 

Figure 2 shows the dependence of the SHG intensity on the applied bias measured at the 
maximum of the rotational anisotropy. Clear oscillatory-like behavior is observed. According 
to Eq.(|I]) this is related to the oscillatory dependence of the DC-electric-field on the bias. 

The following model is developed to obtain the DC-field distribution. First, we suppose 
that there is no charge coupling between the upper Si layer and the metal electrode. There- 
fore, the current through the structure is zero, the Fermi level is the same for all QW's and 
the Si substrate. Second, since {ep/ec)± << 10~^ the structure should be considered as a 
set of charged QW's, with the 2D Fermi-gas of electrons in each well. 

We consider the dependence of E^'~' and the electrostatic potential (f in the middle of 
the Si layer and the 2D charge density n as functions of the discrete variable the number 
of well /. The following discrete analogue of the electrostatic equations with appropriate 
boundary conditions is used: 

' -(E^^(Z) - - 1)) = -27r(n(/) + n{l - l))/e 

< - 1) - = (ed + + E^^{1 - l))/2e , (2) 

ip{0) = 0,v{L) = Uo 

where d and D are the thicknesses of the Si and Si02 layers, respectively; e and ei are static 
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dielectric constants of these layers, Uq is a bias voltage. 

We take into account only the first quantized sub-band of the electron energy spectrum 
of the quantum well. The expression for n{ip) is obtained by the same procedure as for the 



3D case [TO and is given by: 



n{ip) = ^—j^ ln[l + exp((ee</9 + fi)/T)] H -3— ln[l + exp{-{eh(p - fi + A)/T)], (3) 

irh Tin 

where /i is the Fermi energy, A is the bandgap, T is temperature, m is the effective mass of 

the carrier, subscripts e and h stand for electrons and holes. 

In the limit of T = and A = the system (2) is linear and ip{l) takes the form: 

^{l) = C,A' + C2A-\ (4) 

where Ci and C2 are constants and 

A=[l- ^{ed + eiD) /e^ao^ (l + ^ {ed + e^D) /e-^a^ ' . (5) 

denotes the Bohr radius Ti^ {me^)^^ . If {ed + eiD) > e^ao, then A is negative and ip{l) 
reveals the oscillatory dependence. This type of dependence remains for the case of the 
finite magnitude of bandgap and temperature, as Eqs. (2) should be solved numerically. 
Thereafter the SHG bias dependence can be obtained straightforwardly. As discussed below, 
the oscillations in ip{l) result in oscillations in E^'-''{Uo) for the buried 2/m layers, which 
give rise to the SHG signal, accordingly in Eq.(l). 

The total SHG intensity is composed of the DC- field-induced dipole and field-independent 
quadrupole terms (Eq. 1). The relative values of the amplitude and phase of the quadrupole 
term as well as the fiat-bend potential are used as adjustable parameters. The best approx- 
imation is shown in Fig. 2 by solid curves with arrows indicating the bias interval where 
only the first quantized electron sub-band is occupied. Outside this interval the potential 
dependence of carrier density n{ip) should be modified with the account of the occupation 
of the second sub-band. The shift of the bias region of the model validity is a consequence 
of the initial bend banding in the silicon substrate. For the thicker structure the adjusted 
position of the fiatband potential is naturally shifted to the higher voltage 0. 
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In discussing the origin of the oscillations, it is useful to compare the system studied 
above with the classical plasma. In the latter the spatial distribution of charge density 
is monotonic, and the screening is well described by the kinetic equation for the electron 
distribution function. This description by the distribution function requires averaging of 
fields over the elementary volume which includes a large number of particles (the so-called 
physically infinitesimal volume). This procedure allows the replacement of the field acting on 
a particle by the average field Jl^ . Generally, the mean- field potential acting on a particular 



carrier is formed by all charges in the system except the carrier itself. This exclusion is not 
important for the 3D Fermi-liquid, in which electron wave functions are delocalized in space. 
Another situation occurs in the MQW structure where the carriers are confined inside the 
QW's. This results in the distinct difference between the acting and average fields. We 
consider the particular case when the Fermi energy for the 2D electron liquid in QW's is 
much smaller than the DC potential. This means that the interaction between electrons of 
the same QW is negligible in comparison with their interaction with other charges of the 
system. Consequently, the acting potential for the electrons in l-th well is the DC-potential 
formed by external charges all the wells, except the l-th. On the other hand, the average 
potential is still contributed by all the wells. 

Figure 3 demonstrates the appearance of oscillations in f{Uo) due to the discussed effect. 
At small bias the occupation of all sub-bands in all QW's is zero (thin solid line). Consider 
the potential Uq^i, when electrons accumulate in well 1. The potential which determines 
the occupation number is purely the DC potential ipei of charges at electrodes from the bias 
supply (the medium dashed line). This is worth noting again that the charges in well 1 do 
not contribute to this potential. On the other hand, they do contribute to the complete 
DC potential in the system (medium solid line). The kink in this line is determined by 
the charge of well 1, i.e. by cpei. For the certain layer thickness this kink can result in the 
negative in well 1. Similarly, at higher potential f/0,2 holes are accumulated in well 2, 
which is illustrated by the solid thick curve. One can note from Figure 3, that ip and i?^*^ 
in buried layers appear to be oscillatory functions of bias Uq in the considered situation. 



In conclusion, the DC-electric-field screening in the Si-Si02 MQW layered structure is 
studied both experimentally and theoretically. The dependence of DC-electric-field E^'^ 
on the applied bias voltage probed in the deep buried Si wells by DC-electric-field-induced 
SHG shows essentially non-monotonic, oscillatory-like behavior. This behavior which is un- 
expected for 3D Fermi liquids, indicates the importance of the correlations in electron sub- 
system for the strongly anisotropic layered system. The observed phenomenon is explained 
using the theoretical description, which takes into account the strong electron Coulomb lo- 
calization in the quantum wells due to the high anisotropy of the effective mass. Because of 
this localization the electron subsystem of MQW differs from the 3D Fermi liquid essentially. 

We are thankful to V.N. Golovkina for the assistance in experiment and P.V. Elyutin 
and A. A. Nikulin for useful discussions. This work was supported by RFBR grants 97- 
02-17919, 97-02-17923, Special grants for Leading Russian Scientific Schools 96-15-96420 
and 96-15-96476, RFBR-DFG grant 98-02-04092, INTAS-93 grant 0370(ext), INTAS grants 
YSF-98-9, YSF-98-10, ISSEP grant d98-701. Programs "Center of Fundamental Optics and 
Spectroscopy" and " Physics of Solid Nanostructures" . 

Figure Captions 

Fig. 1. Upper panel: the schematic of MQW structure with vicinal substrate, and the 
configuration of the experiment. Lower panel: Azimuthal dependences of the SHG intensity 
for MQW with d — 1.1 nm and D — Bnm, measured in transmission geometry 
for parallel polarizations of the fundamental and SH waves, for two biases: (•) - -1-2.6 V, 
(o) - +1.26 V. Solid lines are fits to data by Eq.(l) with components of x*"^*^"* cind x^^^\ 
corresponding to the 2/m symmetry, as adjustable parameters. Upper panel: schematics of 
the sample and the geometry of experiment. 

Fig. 2. Bias dependences of the SHG intensity for MQW with D = 3nm (upper panel) 
and D — Bnm (lower panel). Solid lines are fits by the model. Arrows indicate the region of 
validity of the model, which is determined by neglecting the second and upper sub-bands. 

Fig. 3. Lower panel: the schematic of the spatial distribution of potential (solid lines) 
across MQW for three monotonically increased values of Uq: C/o,o, t^o,i, t^o,2- The second 
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subscript indicates the number of wells participating in the screening of the external DC- 
electric field. Directions and relative values of DC-field in buried layers corresponding 
to these values of Uq are shown by arrows. Dashed lines are ipei (see text). The upper panel 
shows the theoretical dependence of the DC-electric-field in buried 2/m layers E^'^ as a 
function of the external bias Uq for MQW with d — 1.1 nm and D — 5nm. 
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